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Abstract 



anticipated luminosity of DA$NE. 



We estimate the experimental accuracies to which the K e 4 param- 
eters will be measured at DA$NE. We consider the accuracies ob- 
tainable using asymmetries, and the maximum likelihood method. 
We find the current determinations of the relevant parameters will 

^ . be improved by a factor of five to ten after a year of running at the 
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1. Estimating experimental uncertainties: Generalities 

In this contribution to the DA$NE physics handbook, we estimate the accuracies 
achievable at DA$NE in measuring the theoretical parameters describing the decay K + — > 
7r + 7r~e + z/ e . We refer the reader to the contribution of Bijnens, Ecker, and Gasser for a 
detailed discussion of the theory of this decay. We should note that we refer to this decay 
as A" e 4; we consider only this particular decay and not the channels with neutral pions or 



kaons. 
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We will consider two of the possible ways of extracting the parameters of a theory 
from a set of data. The first one is the classic technique of asymmetries. K e 4 decays have 
a partial decay rate d 5 r of the form 

d 5 T = G 2 F \V us \ 2 N(s n ,s l )J 5 (s n ,s l ,9 7V ,9 l ,4>)ds 7 ,ds l d(cos9 7V )d(cos9 l )d(t). (1.1) 

Here s^ and si are the effective mass squareds of the dipion and dilepton systems; 9^ is 
the angle of the it + in the dipion cm. frame, with respect to the dipion direction in the 
K cm. frame, Q\ is the analogous angle for the e + , and <fi is the angle between the pion 
plane and the lepton plane in the K cm. frame. A is a kinematical factor depending 
only on s^ and S[ and none of the theoretical parameters. The quantity J5 can be written 
as an expansion in sines and cosines of and 9[ multiplying nine intensities I{, which can 
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in turn be written in terms of the three remaining kinematical variables, and three form 
factors, F, G, and H. The explicit dependences may be found in Ref. 1. As can be seen 
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in the contribution of Colangelo, Knecht and Stern, the tangent of the phase shift, 
tan(5o — Si), can be neatly extracted from the ratio of the intensity functions, J7/2/4, or, 
equivalently, 17/2/4, where by the tilde we denote intensities integrated over s n and si as 
well as over 9^. All of the ij, in turn, can be written as asymmetries; in particular we 
have 



T« h 



( <k/2 5tt/2 


tt/2 3tt/2 
r r 


ty 5tt/2 tt 3tt/2\ 




r 


//- 


//" 


7/ + // 


sin^d^d^ 


J^ds^dsidO^ 


\0 3tt/2 


J J 

tt/2 


7r/2 37r/2 tt/2 tt/2 / 


J 





1.2) 




and 

I7 — I / ~ / I d(f) J J5ds 7r ds/dcos^ 7r dcos6' Tr . (1.3) 

\0 7T / 

The asymmetry presents a transparent, elegant and quick (computationally) way to de- 
termine a parameter. 

The second, and main, method we consider is that of the maximum likelihood, which 
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we shall refer to as the MLM . Let x be the vector of phase space variables specifying 
an event in an experiment, and f(x;p) be the probability distribution function predicted 
by a theory, p is the set of parameters in the theory, to be determined experimentally. 
The probability of observing an event at x in the interval d n x is f(x;p)d n x. The function 
/ is normalized to 1 over the whole x interval in which x is physical. In particular, if cuts 
are imposed on the phase space, / must be normalized over the reduced phase space; we 
must also be certain that the normalization is maintained even when the parameters are 
varied from their central value. 

The likelihood, or joint probability distribution, of an experiment yielding iV events, 
each specified by a set of phase space variables Xi , is then defined as 

jV 



C = Jlf(x i ;p). (1.4) 
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The best estimate for p is then simply the value p which maximizes C (or equivalently, of 
W = log£, which is easier to compute). It can be shown that the error matrix, in general 
non-diagonal for correlated parameters, is 



(Pi~Pi)(Pj-Pj) = -^ 



1 fdf(x;p)df{x;p) , 



/ V dpi dpj 
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Here the [ j 1 denote matrix inversion. This integral is, in general, not possible to com- 
pute analytically, but is easily evaluated numerically. If we start with a non-normalized 



probability function V, then the following equation is useful: 



1 df df m 1 da da [I dV dV m 

-d n x = — = — h / -^^— -t— d n x, (1.6) 



/ dpi dpj a 2 dpi dpj J aV dpi dpj 

where a = J Vd n x. 

While this method is more complicated and time-consuming computationally, it has 
the decided advantage that it yields the absolute best possible determination of any given 
set of parameters. In today's age of fast computers, it is thus the method of choice. It 
also allows one to determine any parameters chosen, while the asymmetry is only good 
for certain parameters. Finally, the experimental uncertainties in the parameters can be 
determined reasonably easily, without need for an actual simulation of the determination 
of the parameters themselves using this method. 

It is simple to see that the MLM should give a better determination of the parameters 
than an asymmetry: while the asymmetry only uses the information of whether an event 
is in one or another half of phase space, the MLM benefits from the information of the 
precise position of the event. For example, consider a process specified by the probability 
distribution 

/(0;a) = (l + acos0)/2. (1.7) 

The parameter a can be determined by the ratio (N\ — A^)/(iVi + N2) = a/2 where N\ 
is the number of events with 9 between and 7r/2, and N2 is the number of events with 
9 between ir/2 and n. For a small compared to 1, the error on a in this determination is 
thus 2/yN where N is the total number of events. For the MLM, the integral in eq. (1.5) 

is 

1 

r x 2 1 

lim / — -dx = - (1.8) 
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and the error on a in this determination is thus y/3/N. While this improvement may 
not seem very impressive, larger improvements may be expected as the parametrization 
becomes more complicated. Using the MLM directly on the parameter that interests 
us, we are sure to use all possible information, and to take into account the effect of all 
possible cancelling uncertainties. We observe an improvement as above, of about fifteen 
percent, in the relative error AJ7/J7 when we go from the asymmetry method to the 
MLM. However, the relative error on S = Sq — Si when S is determined directly via the 
MLM is two-thirds of its error when determined indirectly by the ratio I7/2I4, if these ij 
are determined using the MLM, and their errors are then combined as uncorrelated. This 
appears to be due to the information lost in integrating the I{ over s^, si, and 9\ before 
using the MLM. 



2. Estimating experimental uncertainties in K e ^ 

For our first estimates, we consider only the first order terms in a partial wave 
expansion of the form factors F, G and H, i.e., we take 

F = f s e iSo G = ge iSl H = he iSl . (2.1) 

This is consistent with the parametrization used by Pais and Treiman, and Rosselet 
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et al. They consider higher order terms, but the coefficients of these terms are found 
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to be consistent with zero by the experiment of Rosselet et al., so we do not consider 
further terms in our initial estimates. 

All the <p and 9i dependence in the problem is contained in the expression of J5 as 
a function of these variables and the If, for this reason these two variables are referred 
to as "trivial." 9^ appears only in the equations for the I{ in terms of F , G, and H (and 
the phase space expansion, should we consider higher order terms). This leaves then f s , 
g and h with a possible dependence on s n and S[. The im phase shifts Sq and 5\ (which 
at this order appear only in the combination 5 = 5o — Si) depend only on s„-. For the 
moment we parametrize f s , g and h by an expression of the form 

y(s n , Sl )=yo(l + \q 2 ) (2.2) 

where q 2 = (s n — Am 2 )/Am 2 , and y stands for /, g, or h. We take the slope A to be the 
same for f s , g and h, and no slope in si at this stage, again consistent with Ref. 5. For 
the dependence of 5 on s^ we will consider average values in a set of 5 bins in s^, and 
consider parametrizations of S in a later section. 

One last important detail remains to be mentioned. The MLM (or asymmetries, 
for that matter) says nothing about an overall factor in the intensity, as we require the 
probability density to be normalized to one. Thus we divide out f S Q (/o for short) from 
the amplitude, as it is the parameter with the most effect on the integrated intensity. 
Wherever go and ho appear, they are divided by /o, so we replace them by new parameters 
<?o = 9o/ fo an d h' = ho/ fo (A and 5 are unaffected). We then apply the MLM to the set 
of parameters S, g' , h and A, and obtain the correlation matrix 
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The diagonal entries of this matrix are variances of the four parameters, where N is the 
number of events. The off-diagonal elements represent correlations between the parame- 
ters; in this case they are small, but they can be significant, depending on the parametriza- 
tion used. We do not report the full correlation matrix for each parametrization in this 
paper, but they are available from our programs if needed for further calculations. They 
can not be neglected in general if one wants to calculate functions of the parameters 



we use, and propagate the errors correctly. In the end of this paper, we consider the 
determination of some highly correlated parameters. 

We then extract the error on fo from the equation 



(2.4) 



rv e4 = en vd n x 



where V is the unnormalized probability function inputted to the MLM calculation, and 
C represents all the constant factors (masses, two's, 7r's) needed to complete the equation. 
The relative error AT/T is given by ps \/y~N in an experiment like KLOE, where the 
statistical error will be dominant. The relative error on the integral (a) is given by the 
matrix product 



Aa / 1 da , . da 
a 2 dpi %J dpj ' 



(2.5) 



Combining this error in quadrature with the statistical error on T, we obtain the error we 
quote for fo; combining the error on fo in quadrature with the errors on g' Q and h' Q , we 
obtain the errors we quote for go and /lo- 
in Table 1 we display the results of this calculation. The central values (our input) 
are those found by the previous experiment. We have used the program VEGAS to do 
the necessary integrals in five-dimensional phase space. The normalization of the proba- 
bility distribution is ensured automatically by the program, and the necessary derivatives 
also computed numerically. Estimated errors are shown for N = 30000 events, the statis- 
tics of the previous experiment) and N = 300000 events, the anticipated statistics in 
one "year"= 10 seconds of running with C = 5 x 10 cm s . All errors in this paper, 
unless otherwise noted, are statistical errors and can be simply scaled by 1/yN for dif- 
ferent numbers of events. As a general rule, also, the fractional error, on parameters such 
as fo, 90 an d ho, is roughly independent of the central value taken for them, while the 
absolute errors on parameters such as 5 and A remain constant. The last line of Table 1 
shows the errors on these parameters found by Rosselet et al.. We do not quote the error 
on 5 because an error on 5 averaged over the whole of phase space is not very meaningful. 
The errors shown are independent of the central value of S used. 



Table 1. Central values and estimated errors for /q, go, Hq and A 





fo 


90 


h 


A 


Central values 


5.59 


4.77 


-2.68 


0.08 


Errors (N = 30000) 


0.029 


0.059 


0.44 


0.011 


Errors (N = 300000) 


0.009 


0.019 


0.14 


0.004 


Errors (Rosselet) 


0.14 


0.27 


0.68 


0.02 



At this point, before going on to further discuss errors in KLOE, it is necessary to say 
some words about why our estimated errors at "Rosselet statistics" are so different from 
those that Rosselet quotes. The errors given above are purely statistical, but apply to a 
"perfect" detector, i.e., one which covers the whole of phase space with unity efficiency 
everywhere. This is close to true for KLOE; we will attempt to illustrate this later in 
this paper, and will describe a more rigorous demonstration in a future paper. However, 
Rosselet's detector was far from "perfect." In the error we have quoted for /o, the errors 
from Tx ei and the normalization a contribute about equally; the first is about 0.6% and 
the second about 0.9%. Rosselet, however, quotes a relative error AIV e4 /IV e4 of 4.5%, 
which completely accounts for their large error on /q. Their fixed target experiment had a 
10% overall efficiency for K e ^, and a highly variable efficiency as well, varying, for example, 
smoothly from > 95% in a very small portion of phase space with large s; and small s n , 
to near zero at large s^ and small s/. KLOE is in contrast a hermetic detector, operating 
at a e + e~ collider running at the resonance, producing self-tagging low momentum K 
pairs. It will have a uniform near-100% efficiency over all of phase space, minus a few 
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percent of phase space that will be cleanly cut and discarded. 

The next step in our analysis was to drop the slope parameter A and determine 
the errors on the parameters in five bins in s % , chosen so as to have equal numbers of 
events. Such an analysis with real data would have the advantages of studying the s n 
dependence in a more parametrization independent way. If, however, the s n dependence is 
correctly given by eq. (2.2), this method will not determine A as accurately, so in general 
both types of approach are necessary. For our purposes, displaying the error in bins is 
also important to illustrate the possible accuracies with which S(s^) may be measured, 
before we implement a possible parametrization of S. In Table 2, we give the estimated 
errors, taking an average of yo(l + Xq ) in each bin as our inputs for y — f, g, and h, 
with Rosselet values for the yo and A. The errors on S are essentially independent of the 
inputs of its central value. In the last line, for comparison, we display the errors on S as 
measured by Rosselet et al. The improvement is not as drastic as that of fo was, but is 
nonetheless a factor of 1.5 to 2. This should be further multiplied by a factor of V^IO to 
V20 per DA$NE running year. The accuracy on / in bins is even better than we might 
have expected from the error on fo multiplied by v5. This is because the error on A 
gives most of the contribution to the error on the normalization a, and thus a significant 
contribution to the error on fo- 

We have examined in some detail the question of what accuracy S can be measured 
to. We have first of all determined that while S appears in I\, I2, I4, I5, I7, and Ig, it 
is only the dependence of I7 that gives us the above accuracy on S. This can be seen by 
replacing S in all the I{, except I7, by a dummy variable So, set equal to the central value 
of S. When we proceed to apply the MLM to the new probability function, we find the 
same error on S as before, within a few percent. If, however, we apply the MLM to the 1% 
as parameters in their own right (we cannot use the Ii as parameters, because they are 
functions of the phase space variables) and then take the ratio to determine S, we find that 



Table 2. Estimated errors in five bins of 6000 events each. 



v^ (GeV) 


0.279-0.3 


0.3-0.316 


0.316-0.334 


0.334-0.357 


0.357-0.494 


/ 


0.037 


0.039 


0.041 


0.043 


0.047 


9 


0.195 


0.142 


0.124 


0.116 


0.112 


h 


1.46 


1.04 


0.93 


0.90 


1.00 


5 


0.062 


0.041 


0.034 


0.029 


0.025 


S (Rosselet) 


0.13 


0.07 


0.05 


0.04 


0.04 



the error on 5 increases by 50%. (If we use the asymmetry method to determine the I{, 
the error increases another 15%.) We have not taken care to cancel correlated errors in I4 
and I7, but we have checked that the correlated parts of the errors are small relative to the 
uncorrelated parts. So, this 50% increase appears to be mainly due to the information lost 
in integrating the I{ over three out of five of the phase space variables before applying the 
MLM. Equivalently, the better error on 5 can be attributed to applying a more detailed 
parametrization (therefore more information) from the beginning of the calculation. 

Nonetheless, it may be interesting to determine the h and their errors as a parame- 
terization independent way to present the data. We have estimated that the combination 
I\ ~ h/S can be determined to 0.4% in five bins of 60,000 events each. The other I{ can 
be determined with absolute errors of one to two times this error. 



3. Other experimental uncertainties 

So far we have estimated the statistical experimental uncertainties in a perfect de- 
tector. In this section we would like to say a few words about other experimental un- 
certainties. Errors (other than those already estimated) in a detector like KLOE, where 
the efficiency is essentially 100% outside of a small region of phase space to be removed 
by cuts, fall into two categories: those associated with cuts, and those associated with 
resolution. We intend to write an event generator, that would allow us to impose a proper 
cut on true kinematical variables such as the angles of the particles themselves, or their 
momenta, but for the moment are only calculating the MLM error estimating integral 
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over si, Sjr, 61, 6 n and 0. We expect that appropriate cuts will be applied to reject 
charged particles that have momenta less than 20 MeV, and those that are within 9° of 
being parallel to the beam pipe (because they will not cross enough wires). We expect 
that both of these cuts will (from the point of view of our problem) reject events fairly 
randomly, and thus only have the effect of reducing the number of events by a couple of 
percent, having a negligible effect on our error estimates. Just to check, we have con- 
sidered cutting s n in the range 4m 2 to 4m 2 + 4 * (20MeV) 2 , and si in the range to 
(20 MeV) 2 . The s n cut has negligible effect and the s/ cut increases our errors by about 
10 percent or less. Eventually we would also like to investigate the effect of inaccurate 
cuts (for example a cut that is believed to be at 20 MeV and is really at 15 MeV), which 



can be an important source of systematic errors. 

Smearing, also known as convolution with a resolution function, can be imposed 
on our variables with more confidence, being a random effect. We assume a gaussian 
resolution function, i.e., we replace the probability function f(x;p) by a new function 
g(x;p) 




9(x;p) = f(x";p) -= dx" « / f(x-x';p) T^dx' 

J o"V 2tt J aV 2tc 

(3.1) 

OO OO 11 v ' 

f(x-x';p)e~ x ' 2/2a2 e- y2/2a2 -^dx'dy = 1 1 f(x - x';p)dadb 

-cx^oo 

where y is a dummy variable introduced to make the gaussian integrable, a = e~( x +y '' 2a 
and tan(27r6) = x'/y. x$ and x\ are the lower and upper bounds on x; x\ — xq is assumed 
to be large compared to a. Thus we have 



showing the equivalence of convolution with a gaussian and integrating after a gaussian 
smearing of the independent variable. Note that the quantity \/—1 log a sin(27r6), for a 
and b uniformly distributed in the interval to 1, is gaussian distributed with variance 1. 

KLOE expects a resolution of about half a percent to one percent in angles and 
momenta. We make the assumption that smearing by one percent (of 2tt or it in the 
case of our angles) in our five variables is a reasonable and generous approximation to 
smearing in the actual kinematic variables. We observe negligible effect on our errors; in 
fact smearing of as much as 5% of the maximum of our variables has no effect except in s^, 
where it yields a ten percent (fractional!) increase in our errors. This is not surprising. 
If one knows one's resolution function, one can compensate for it, and the accuracy is 
unaffected, at least if the scale of the smearing is small compared to the scale of the 
effect (this scale in our case is of the order of the whole of phase space). If however 
there are unknown parts to the resolution function, i.e., systematic errors, these can result 
in systematic errors in the result, that the MLM error estimating technique above will 
never find, simply because we have no way to input an unknown error. These errors 
could in principle be estimated by applying the MLM itself to simulated data, but this 
is beyond the scope of our investigation. The accurate estimation of systematic errors is 
at any rate something that will have to be done by KLOE. Meanwhile, however, we have 
made one attempt to estimate the effect of such systematic errors by returning to the 
determination of 5 from the asymmetry method. Here we can easily calculate what 5 we 
would "measure" from an imaginary set of data described by a set of input parameters; 
normally of course, we get back the value of S we input. If, however, we replace the 
probability distribution f(x) by a convoluted distribution g(x), then, by the shift in the 



recovered 8, we can see what would be the effect of a resolution function that we did not 
know about, and therefore did not compensate for. 

We have first verified that gaussian smearing still has no significant effect, which 
means that it is not neccessary to know the exact form of the resolution function, as long 
as it is symmetric and not too broad in the appropriate scale. We have next examined the 
effect of an actual shift in each of our phase space variables. The effect is tiny compared 
to the statistical errors for shifts in $i, 6 n , s/ and s n . For a 1% shift in 0, S goes down 
by 0.01, or about half of the statistical error with 300,000 K e 4 events. We remark that a 
1% systematic shift in <fi is already so large as to be inconcievable, as is a difference of 
observable angles. 

4. Other parametrizations 

In this section we would like to present results for some extensions of the minimal 
parametrization we have used so far: (1) we include a p-wave term in F, and simulta- 
neously adopt a notation consistent with the partial wave expansion of Ref. 1; (2) we 
consider an si dependence in /, g, etc.; (3) we consider two parametrizations of S in terms 

l Q 

on this parameter. 

F thus becomes 



of s-wave scattering length aS to give a first indication on what errors can be expected 



F = f s e iSo + f p e^ cos9 n - ^lL± cos9 n ge iS \ (4.1) 

where a n , P ■ L and X are defined in Ref. 1. Our previous parametrization is thus 

equivalent to this one for f p = <J 7T P- Lg/X. For the s/, dependence, chiral perturbation 
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theory to one loop predicts the same slope for s n and s/, or, more precisely, 

y(**>'i) = yo(i + *( Sw ~t™l~ sl ))- (4-2) 

To be more model-independent, and to have a separate determination of the s^ and si 
slopes, we choose the parametrization 

y(sn, si) = y (l + \q 2 ) (1 - \2Si/4ml) . (4.3) 

We then find, for 300000 events, 

A/ s0 = 0.014 Ag = 0.038 Ah = 0.14 Af p0 = 0.014 (4.4) 

AA = 0.004 AA 2 = 0.011 (4.5) 

to be compared with table 1. We have displayed the errors for the old parameters as 
well as the added ones, as some of them have increased slightly. We have taken central 
values A2 = 0.08 and f p o = 3.3 (which maintains the normalization unchanged); if other 
parameters are preferred it should be remembered that the fractional error of f p o, and 
the absolute error of the slope, remain essentially constant. 



For S we have first considered the parametrization used by the previous experi- 
ment, to compare our estimated errors with the ones they determined. They use 




sin2 5 = 2 J Jn 4m " (og + bq 2 ) , (4.6) 

y s ty 

where b = 6q — a\, the difference between the s-wave slope and the p-wave scattering 
length. There is moreover a possible relation between b and a[] 

6 = 0.19-(a[]-0.15) 2 . (4.7) 

For 30000 events (to compare with Ref. 5) we find Aoq = 0.029, to be compared with 0.05 
in Ref. 5, if we use both eqs. (4.6) and (4.7). If we use only eq. (4.6) we find Aajj = 0.06 
and Ab = 0.07, to be compared with 0.11 and 0.16 in Ref. 5. We have used the central 
values aj] = 0.28 in the first case, and a[j = 0.31, b = 0.11 in the second case, as found in 
Ref. 5. So, in one year running at DA$NE at a luminosity of 5 x 10 32 cm 2 s _1 , we expect 
a factor of five improvement in the error on the mr scattering length, meaning that the 
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existing discrepancy between measurements and predictions for aR should be resolved. 



We have also considered a more recent parametrization, due to Schenk. He gives 



tan tf (*„) = J^i^ (a{ + b\q 2 + c\A ( ^L^L ) (4 . 8 ) 

where q 2 = (s n — 4m 2 ) /Am 2 and 

~ r r t AIT! / r \ 3 

Here / denotes isospin, I angular momentum, and 5{j is the Kronecker delta. We need 
to calculate 6q = 5® and 5\ = 8\\ Sq and s\ (the values at which the phase shifts should 
pass through 90°) are the squares of the a and p meson masses. There are far too many 
parameters here to be able to determine them by K e ^ measurements alone; to begin with, 
we fix the higher order coefficients c[] and c\ to zero. The remaining set a§, a\, 6q and b\ is 
still too large: if we plot 5 versus s n , we find that while changing a\ from its central value 
of 0.037 to a new value of 0.087 changes S essentially uniformly by 15 percent, changing 
instead 6q from its central value of 0.24 to 0.19 has the same effect, to within about 1 
percent. Moreover, with appropriate variations in a[] and b\, one can quite well mock up a 
variation in a\ (or &[], naturally). Thus, one can only hope to determine two independent 
parameters without recourse to other experiments or other theoretical constraints. 

The MLM quantifies these conclusions. Assuming Schenk's central values (oq = 0.2, 
b\ = 0.005, and above), Oq or b\ can be determined to 0.03 (for Rosselet statistics of 30,000 
events), or 6q or a\ to 0.04, if only one parameter is considered, and the rest held fixed. 
If one determines a[] and b\ simultaneously, holding the rest of the parameters fixed, the 
correlation is 0.71 (defined as the covariance of the two parameters divided by the cr's 
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for each, 1 for a maximally correlated pair of variables), and the errors are 0.04 and 0.05 
respectively. If one determines instead a[] and 6q — a\ simultaneously, the correlation is 
—0.9 and the errors are 0.06 and 0.08; for b\ and 6q — a\ the correlation is 0.94 and the 
errors are 0.09 and 0.11. To determine more than two parameters simultaneously is not 
possible. 

5. CONCLUSION 

We have considered asymmetries and the maximum likelihood method, and shown 
how the latter is more powerful and yields more precise determinations. We have found 
that the current measurements of these parameters will be improved by a factor of five to 
ten in one effective year (10 7 seconds) of running with C = 5 x 10 32 cm 2 s _1 ; in particular 
we expect that the irir scattering length a§ should be measureable to an accuracy of about 
0.01. 
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